Removal of interference from gravitational wave spectrum 
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Abstract. We develop a procedure to remove interference from gravitational wave 
spectrum. The method is applied to the data produced by the Glasgow laser inter- 
ferometer in 1996 and all the lines corresponding to the interference with the main 
supply are removed. 

1 Introduction 

In this paper we present an algorithm to remove interference from the gravita- 
tional wave (gw) spectrum. This method allows the removal of coherent lines 
coming from deterministic signals while keeping the intrinsic detector noise. 
Unlike other existing methods for removing single interference lines ^ |tJ, the 
method described here can remove the external interference without removing 
any gravitational wave signal that may be hidden by the interference. There- 
fore, it can be very useful in the search for monochromatic GW signals as those 
ones produced by pulsars H || . 

The key to this method is to determine the interference by using many har- 
monics of the interference signal. In the study of the data produced by the 
Glasgow laser interferometer in March 1996 ||, we observe in the spectrum 
many instrumental lines, some of them at multiples of 50 Hz. All these lines 
are wide, and when we compare them, we observe that their overall structure 
is very similar but only the scaling of the width is different (see figure 1). If 
we look at these lines in more detail, in smaller length Fourier transform (sec- 
onds in length) , they appear as well defined small bandwidth lines which change 
frequency over time in the same way, while other ones remain at constant fre- 
quency. Therefore, all these lines at multiples of 50 Hz must be harmonics of a 
single source, for example the main supply. 

In the Glasgow data, the lines at 1 KHz have a width of 5 Hz. Therefore, 
we can ignore these sections of the power spectrum or we can try to remove this 
interference in order to be able to detect GW signals hidden behind those lines. 
The LIGO group has also reported largely instrumental artifacts at multiples of 
the 60 Hz line frequency in their 40 m interferometer jlj. Therefore, this seems 
to be a general effect present in the different interferometers. 

This electrical interference may possibly be reduced by improvements in the 
electronic design that will be incorporated in the next generation of detectors. 
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Figure 1: Comparison of the structure of the lines at 250 Hz and at 350 Hz 
of the power spectrum of the Glasgow data. The broad shape is due to the 
wandering of the incoming electricity frequency. 

However, we cannot be sure that such interference will be completely absent or 
that other sources of interference will not manifest themselves in long-duration 
spectra. Indeed, the Glasgow data || contain other regular features of unknown 
origin Q. For this reason we investigate solutions to the problem using the 
existing data. 

We address here only the removal of coherent lines, not stochastic ones 
(such as those due to thermal noise). Our method requires coherence between 
the fundamental and several harmonics. If there is no such coherence, other 
methods [|[ [?J can be used, but these will remove GW signals too. The method 
presented here is the only one we know that allows one to find real signals under 
the interference. 

2 The basic concept 

We suppose that a electrical signal at ~ 50 Hz enters in the electronics of the 
interferometer. The instrument then transforms the input signal into a series of 
harmonics in a stationary way. Therefore, the data stream contains the function 
y(t) given by 

y(t) = a n m(t) n + (a n m(t) n )* , (1) 



where a n are constants and m(t) is nearly monochromatic around 50 Hz. The 
total output, x(t), contains also the gravitational wave signals s(t) and the 
detector noise n(t) 

x(t) = y(t) + n(t) + s(t) . (2) 

We know that the data recorded is band-limited since an anti-aliasing filter 
was applied to the data before it was sampled. Therefore, the function y(t) 
must be also band-limited, and hence, the sum can be limited to the first 39 
harmonics. This number is given by the Nyquist frequency that is related to 
the sampling frequency of the experiment (n max = I Nyquist/ 50 Hz — 1). 

Our purpose is to construct a function 

39 

Kt) = PnM{t) n + ( Pn M(t) n )* , (3) 

71=1 

similar to y(t). Thus, we have to determine the complex function M(t) and all 
the parameters p n . Notice that from the experimental data we do not indepen- 
dently know the value of the input signal m(t). 



3 The algorithm 

In this section we present the algorithm to remove interference from the gravi- 
tational wave spectrum. 

As we pointed out in the previous section, we assume that the data produced 
by the interferometer is just the sum of the interference plus the detector noise 
and possible GW signals that we will not consider here 

x(t) = y(t) + n(t) , (4) 

where y(t) is given by Eq. ([!]). Therefore, the Fourier transform of the data 
x(v) is simply given by 

x{v) = y(v) +n(u) . (5) 

First, we select a set of harmonics (in our case the odd harmonics were 
strong, i.e., k — [3,5,7,9,...]), or all of them if we prefer. The idea is to 
construct the function M(t) getting the maximum information inbeded in the 
harmonics considered. For each of these harmonics, we determine the initial 
and final frequency of the line (i/jfe, Vfk) and we define a set of functions z,k{v) 
in the frequency domain as 

z ^) = { elsewhere! (6) 
Comparing Eq. (0) with (Q) and (||) we have 



Z k (y) = a k m k + n k {v) , 



(7) 



where 



j n(v) v ik < v < v fk 
n ^) = \ o elsewhere (8) 



is considered to be a zero-mean random complex noise. 
Then, we calculate their inverse Fourier transform 

z k (t) = a k m(t) k + n k {t) , (9) 

and we define 

B k {t) = [z k (t)f k = (a k y/ k m(t)f3 k (t) , (10) 

where 

n k (t) 



l/fc 



(11) 



a k m(t) k _ 

All these function {B k (t)} are almost monochromatic at ~ 50 Hz, but they have 
different amplitude and phase. In order to compare them, we construct another 
set of functions 

b k (t) = am(t)0 k (t) , (12) 

by performing the operation 

b k (t) = T k B k (t) , r fc = ]T SiO'At)B fe (jAi)*/ ^ \B k (jAt)\ 2 . (13) 



These new functions {b k (t)} form a set of random variables —functions of time— 
and they all have the same mean value 

{bi) = am(t) . (14) 

Now, we want to construct M(t) as a function of all {b k (t)} 7 in such a way 
that it has the same mean and minimum variance. If we assume the function 
M(t) to be linear with {b k (t)} 7 the statistically the best is 

"<« = (?VC&) /(?«)' < 15 » 

where 

Varfe(i)1 = .l MM +cort ec ti o„ s , (10) 

a k m(t) k w « fc (t) , (17) 

(n fc (t)n fc (t)*) = ( dv ( dV '(n k {v)h k {i/y)e™ { r- v "> t . (18) 




Figure 2: Decimal logarithm of the periodogram of 128 blocks (approximately 2 
minutes) of the Glasgow data, (a) One of the harmonics near 754 Hz. (b) The 
same data after the removal of the interference as described in the text, (c) The 
same experimental data with an artificial signal added at 452 Hz. (d) The data 
in (c) after the removal of the interference, revealing that the signal remains 
detectable. Its amplitude is hardly changed by removing the interference. 



If we assume the noise to be stationary (i.e., (n(v)n(v')*) = S{v)S{v — z/)), the 
previous equation becomes 



(n k (t)n k (ty) = du dv 1 S k {p)5{v - I /) e *ri("-^)t = / S {v)dv , (19) 



where S[y) is the power spectral density of the noise. 

Finally, it only remains to determine the parameters p n in Eq. (||) that we 
obtain by performing the operation 



We have applied this method to the data taken from the Glasgow laser inter- 
ferometer in 1996 and we have succeeded in removing the electrical interference. 
The same method has also been applied to the true experimental data with an 
external simulated signal at 452 Hz, that remains hidden due to its weakness, 
and we have succeeded in removing the electrical interference while keeping the 
signal present in the data, obtaining a clear outstanding peak over the noise 
level (see figure 2). This will be described in detail elsewhere ||. 

This method can also be applied to any other kind of interference and it may 
have more applications, not only for the detection of GW radiation, but also, for 
example, in radioastronomy. 
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